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Abstract. We discuss the bihamiltonian geometry of the Toda lattice (peri- 
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odic and open). Using some recent results on the separation of variables for 
bihamiltonian manifold, we show that these systems can be explicitly integrated 
via the classical Hamilton- Jacobi method in the so-called Darboux-Nijenhuis 
coordinates. 



1 Introduction 

X 



In this paper we reexamine the classical A„-Toda systems with the aim of show- 
ing that these lattices fall into a notable class of bihamiltonian integrable sys- 
tems: those for which a distinguished set of coordinates (the so-called Darboux- 
Nijenhuis (DN) coordinates) allows the solution of the Hamilton- Jacobi equations 
associated with the Hamiltonian flows by means of an (additive) separation of 
variables (SoV) fllB" , p]]. In particular, we will show that such coordinates arise 



from the geometry of the Poisson pencil after a Hamiltonian reduction process 
on suitable symplectic leaves. 

DN coordinates (see, e.g., [jllj) can be naturally defined on a Poisson-Nijenhuis 
(PN) manifold ||, that is, on a 2n-dimensional manifold M. endowed with a sym- 
plectic two-form <jj and a (1, 1) torsion free tensor N satisfying certain compatibil- 
ity conditions. In M we present and discuss an intrinsic condition to characterize 
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those Hamiltonian systems on A4 for which DN coordinates separate the corre- 
sponding Hamilton-Jacobi equations. Moreover, some of the connections between 
Hamiltonian hierarchies which satisfy a certain recursion property with respect to 
the tensor N, and Gel'fand-Zakharevich (GZ) systems 0] are investigated there. 
This paper is devoted to frame the Toda lattices into such a scheme. 

The plan is as follows: in Section |] we sketch the main points of the above- 
mentioned SoV theory for bihamiltonian manifold, referring to for complete 
proofs and a more detailed discussion. Section |3] contains a formulation of the 
Toda lattice within the GZ scheme, that is, taking as starting point its Poisson 
pencil and the problem of finding the Casimir functions. Section ^ concerns the 
application of the bihamiltonian SoV theory to this family of integrable systems. 
Finally, in Section [| we treat the three-particle case to give a feeling of how the 
method works. 



2 Separation of variables on PN manifolds 

Let (Ai,uj) be a 2n-dimensional symplectic manifold endowed with a Nijenhuis 
tensor field N compatible with uj (in the sense of the theory of bihamiltonian 
manifolds). These manifolds are called Poisson-Nijenhuis manifolds ||. Exam- 
ples of such manifolds are provided by bihamiltonian manifolds endowed with a 
pair of Poisson bivectors (Pq, Pi) one of which, say Po, is invertible, In this case, 
uj = P^ 1 and N = PiP Q l . 

Definition 2.1 By Darboux-Nijenhuis coordinates on M. we mean any system 
of local coordinates (Xj, fij)j=i , n which enjoy the following two properties: 
i) uj takes the canonical form 



uj = 

i=i 



^ dXi A dfii ; 



ii) the adjoint Nijenhuis operator N* takes the diagonal form 
N*dXj = XjdXj , N*dfij = Xjdfij . 



It has been shown [1 1| that DN coordinates exist on any PN manifold where N 
has n functionally independent eigenvalues. In this case the coordinates Xj can 
be computed algebraically as the roots of the minimal polynomial of N, 

C(X)=Det(N-Xl)K (2.1) 

On the contrary, the complementary coordinates must be computed (in general) 
by a method involving quadratures. 

In [Q we characterize a class of Hamiltonians on M. whose associated Hamilton- 
Jacobi equations can be solved by separation of variables in DN coordinates. Let 
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(Hi, . . . ,H n ) be a set of functionally independent (Hamiltonian) functions that 
are in involution with respect to the canonical Poisson bracket defined by 

{f,g} = u(X f ,X g ) . 

We assume that the Lagrangian foliation defined by the functions Hi is invariant 
with respect to N. This is tantamount to saying that, at a generic point mGM, 
the differentials dH^ span an n-dimensional vector subspace of T^M. which is 
closed under the action of N*. Hence there exists an n x n matrix F, whose 
entries are functions on A4, such that 

n 

N*dHi = H dH 3 • (2-2) 

Definition 2.2 The Hamiltonians (Hi, . . . , H n ) are separable in the DN coordi- 
nates if there exists an n x n invertible matrix T and an n-component vector V 
such that 

TH — V , (2.3) 

where H = (Hi, . . . , H n ) T , and the matrix T and the vector V possess the Stackel 
properties: 

1. the entries of the j th row of T depend only on the conjugated coordinates 
(Xj,fij). 

2. the j th component of the vector V depends only on (\j,fij) as well. 
A remarkable "separability test" is given by the following 

Theorem 2.3 The Hamiltonians (Hi, . . . ,H n ) are separable if and only if the 
matrix F verifies the equation 

N*dF = FdF . (2.4) 

Two remarks are in order to explain this theorem. First of all, equation (|2.4| ) 
must be read as follows. In the left hand-side, dF is the matrix whose entries 
are the differentials of the entries of F, and iV* acts separately on each entry. 
Secondly, one should notice that ( |2.4p is a coordinate free test of separability, 
that can be checked without computing the DN coordinates. Once the test is 
passed one can construct the Stackel matrix T, still in general coordinates, by a 
simple algebraic procedure. One has to consider the eigenvectors of the matrix F 
and form with them a (suitably normalized) matrix T that diagonalizes F: 

F = T" 1 AT, A= diag(Ai,... ,A n ) . 
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By condition (|2.4 ), this matrix is a Stackel matrix; by condition fl2.2|) , the vector 



= TH verifies the Stackel property. Then, once constructed the DN coordi- 
nates, the Hamilton- Jacobi equations associated with (Hi . . . , H n ) can be easily 
solved by separation of variables. Notice that the DN coordinates separate at 
once the HJ equations associated with any of the Hamiltonians H^. 

To complete the construction of the DN coordinates, that is, to construct 
algebraically the coordinates /ij conjugated to the eigenvalues Xj of N, the fol- 
lowing procedure is often useful. We consider the Hamiltonian vector field Y 
associated (by the symplectic form uj) with the function ^Tr(N), and the space 
of functions F(x; A), depending smoothly on x G Ai and holomorphically on the 
parameter A. We denote with F(x; Xj) the evaluation of F(x; A) at A = A.,-. If 
N*dF(x; Xj) = XjdF(x; Xj) for all j = 1, . . . ,n, we say that F(x; X) is an exact 
eigenvector of N*. 

Theorem 2.4 If F(x;X) is an exact eigenvector of N* , satisfying the "normal- 
ization property" Y(F(x; A)) = 1, then the evaluation of F(x; X) at the points 
X = Xj, i.e., 

Uj = F(x; Xj) , 
provides a set of n remaining DN coordinates. 

In the application of Section ^, we will use the property that if F(x; X) is an 
exact eigenvector, then Y(F(x; A)) is an exact eigenvector as well. Since in the 
separable suitable combination of the Hamiltonians is exact, one can act 

with Y on such a combination and generate a space of exact eigenvectors where 
the equation Y(F(x; A)) = 1 may be solved algebraically. 



2.1 DN separable Hamiltonians from GZ systems 

Let M. a (2n + /c)-dimensional manifold endowed with a pencil P\ = P\ — XPq of 
Poisson tensors. We suppose that it admits k polynomial Casimir functions 

#(«> = ^H ( f ) X na ~ j , a = l,...,k, 

3=0 

(i) 

with n = rii + • • • + nk- If the functions ffj are functionally independent, then 
Ai is called a complete GZ manifold, and the pencil Pa is said to be a pure 
Kronecker pencil of type {2ni + 1, . . . ,2n& + 1}. Since the functions Hq form 
a maximal set of independent Casimirs of Pq, the generic symplectic leaf S of Pq 
is the 2n-dimensional submanifold given by Hq = C a , for a = 1, . . . , k. The 
restrictions H^ to S, for j a — 1, . . . , n a: and a — 1, . . . , k 7 of the n remaining 
Hamiltonians define a completely integrable system in the Liouville sense. 

In order to solve by SoV this system, we suppose that there exist k vector 
fields Z a (to be called transversal vector fields) spanning a fc-dimensional inte- 
grable distribution Z and satisfying: 
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a) The normalized transversality condition: Liez a (H^) = 8 h a for all a,b = 
1, . . . , k; 

b) The deformation condition for the Lie derivatives: Liez a (P\) = J2b=i ^b A Y^ 
for some vector fields Y b a ; 

c) The "flatness" condition: Lie Za (Lie Zf) (#( c )(A))) = 0, Va, 6, c. 

Conditions a) and b) imply that the distribution Z is transversal to the sym- 
plectic leaves of Pq, and that the functions vanishing along Z are a Poisson 
subalgebra with respect to the Poisson pencil P\. Then, as a consequence of the 
Marsden-Ratiu theorem we have that: 



Proposition 2.5 The Poisson pencil on Ai can be projected on the generic sym- 
plectic leaf S of P , so that S becomes a PN manifold. The functions \ for 
j a = 1, • • • , n a , and a = 1, . . . , k, satisfy the condition ( $.2j , that is, there exists 
an n x n matrix F such that N*dH = FH , where H is a column vector collecting 
the above functions. 



Under the "flatness" condition c), one can show that equation ( |2.4j) is satisfied, 
so that the reduced Hamiltonian system is separable in the DN coordinates. 
These coordinates may be computed from the geometry of Pa, without actually 
performing the reduction process. In this case, in fact: 

1. The minimal polynomial of the Nijenhuis tensor N induced, according to 
the previous proposition, on the leaf S is the determinant of the matrix 
G(X) = [Ue Za (H^(X))] ab=1 k , that is, detG(A) = iff A = A,; 



2. The vector field Y of Theorem |2.4| is given by Y = Yla=i ^o"! 



3. If (1, p 2 (A), • • • , Pfc(A)) satisfies 

(l, P2 (A),...,p fc (A))G(A)=0 forA = A 



■3i 



then #W(A) + p 2 (\)H { - 2 \\) H h p k {\)H^ k \\) is an exact eigenvector 

of iV*. Hence it can be used to find a normalized exact eigenvector, and 
therefore the pj coordinates. 



Remark 2.6 The SoV theory for PN manifolds outlined above provides intrinsic 
and algorithmic recipes to check whether a given Liouville integrable system 
defined on a PN manifold can be separated in the DN coordinates. On the 
other hand, the conditions under which one obtains separable Hamiltonians from 
a GZ manifold are by no means algorithmic. In particular, the existence of 
the distribution Z (that is, of the vector fields Z a and Y% fulfilling the above 
three properties) must be checked (and guessed) case by case. In ||, [3J] some 
GZ systems, obtained from stationary reductions of the Boussinesq and KdV 
hierarchies, are discussed along these lines. In the next sections we will apply the 
scheme herewith outlined to the Toda lattices. 
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3 The Bihamiltonian approach to Toda lattices 



The phase space of the (complex, periodic) Toda lattice (see, e.g., [|]) with n 
sites (particles) is the manifold Ai = (C*) n x C n parametrized by the Flaschka 
coordinates {a,, 6j) i= i,... in . We endow it with the Poisson pencil Pa defined as 
follows (see, e.g., [JTJJ] and references cited therein). It associates with the one- 
form J2 k (akddk + Pkdbk) the vector field ^2 k (a>kd ak + bkdb k ) according to the rule 



dk = a k ((b k - X)(3 k - (b k+1 - X)Pk+i + a k -ia k -i - a k+1 a k+1 ) 



b k = (b k - X)(a k - 1 a k ^ 1 - a k a k ) + a k f3 k+1 - a k _ x ^ t 



(3.1) 



fc-i 



where the cyclicity condition (-) k + n = (-) k is implicitly assumed. We write the 
matrix expression of P\ = P\ — XPq in the 3-particle case, the n-particle case 
being easily generalized from this example: 







-a\Q2 CL1CL3 0.1(61 — A) —0-1(62 — A) 







—0203 



-a 3 (6i - A) 



02(62 - A) -02(63 - A) 
03(63 - A) 



oi 



-o 3 



02 



(3.2) 





According to the GZ scheme, we study the kernel of Pa- We have to solve the 
equations 

(6 fc - X)f3 k - (6 fc+ i - X)f3 k+ i + a fe _ia fc _i - a k+1 a k+1 = 

(6fc - A)(a fc _ia fe _i - a k a k ) + a k (3 k+1 - a fc _i/5 fc _i = 

With algebraic manipulations (see |I3)|), it can be traded for the system of equa- 
tions 



{b k ~ X)/3 k + o fc _ia fc _i + a k a k = L x 
{a k a k ) 2 + a k p k p k +i - L x a k = L 2 



(3.3) 



where Li are Z n -invariant functions. Setting L\ = 1, L2 = 0, and introducing the 
variables 

ilk — , 

"A; 



we obtain the following Riccati type equation: 

hkhk+i = (frfc+i — A)/ifc + a k . 



(3.4) 



Proposition 3.1 The characteristic equation ( \3.4\ ) admits a solution hk which is 
a Laurent series in the parameter A of the form hk = A + ^2JLi ^fcjA~ J . The Lau- 
rent coefficients hkj can be computed by recurrence as functions of the variables 
{<2j, bi}. The product C = hi • • • h n of the components of any solution of (\3.J\ ) is 
a Casimir function of the Poisson pencil P\ . 

Notice that, once the characteristic equation is solved, the one-forms in the kernel 
of Pa can be easily computed (by recurrence) solving the system 



h k a k + a k Pk = 1 

dkhk = Ac+l 



k = 1, . . . , n, 



which is equivalent to the system ( |3.3| ) with L 1 — 1, L 2 — 0. 

This method allows us to find Casimirs of Pa that are Laurent series in A. 
According to the GZ scheme 0, however, we should better look for polynomial 
Casimirs of P\. They can be found linearizing the Riccati equation (|3.4j ) as 
follows. 

Setting hk = fiipk/ipk-i, we transform equation ( |3.4|) into the linear system 



0. 



(3.5) 



where fi is related to the Casimir C via C = [i n . In matrix form we have Lip = 0, 
where 



C 



nibi — A) — fi 2 
cii fi(b 2 - A) 
a 2 













/i(6 n _i - A) -/i 2 

O-n-l M fo n - A) 



(3.6) 



This is how the classical Lax matrix of the Toda lattice can be introduced into 
the game in the GZ bihamiltonian point of view. We remark that, since the 
Riccati equation (|3.4|) admits solutions, so does the linear system Cip = 0. So, 
taking into account the cyclicity of C, we arrive at 



7 



Proposition 3.2 The spectral curve of the problem, det(£) = 0, is a quadratic 
polynomial in the Casimir C , 

det(£) = -C 2 + P (1) (A)C + P (2) . (3.7) 

Thus, both ffW(A) and are polynomial Casimirs of P\. In particular, 

# (2) = (-l) n+ V--a n 

is a common Casimir of Pi and P , and H^'(X) has the form 

n 

H {1 \\) = (-l) n A n + ^Hf ] \ n ' ] . 

3=1 

It can be easily realized that if . = ( — l)- 7 cr"(6 1 , . . . ,6 ra )+lower order terms in 
the bj, where aj is the j—th elementary symmetric polynomial in n letters. So, 

the Hamiltonian functions {H^ 2 \ {if,- } 3 -=i,... )n ) are functionally independent and 
the previous proposition provides another proof of the fact M that the periodic 
Toda lattice with n particles is a complete GZ manifold of type {1, In — 1}. 

We end this section with a remark which frames the open Toda lattice within 
this scheme. It is well known that the open Toda lattice can be obtained form 
the periodic one by pulling one particle to infinity, that is, in the Flaschka coor- 
dinates, by letting one of the a coordinates, say a n , attain the value. The 
phase space of the (complex) open Toda lattice with n particles is thus the 
manifold Ai = (C*) n_1 x C n parametrized by reduced Flaschka coordinates 
{ai . . . , a n -i,b\, . . . ,b n }. The Poisson pencil Pa of the open case can be ob- 
tained from the periodic one by means of the following trick. Let Ai D Ai be 
the manifold obtained from the phase space of the periodic lattice adjoining the 
"boundary component" defined by a n = 0. The Poisson pencil defined by (|3.1|) , 
being polynomial in the Flaschka coordinates, extends naturally to a Poisson 
pencil P\ on the extended manifold Ai. The phase space Ai of the open case can 
be identified with the zero set of the common Casimir of Pi and Pq, which, 
obviously enough, is still given by = (— l) n+1 a! • ■ ■ a n . Then Pa can be re- 
stricted to Ai, and its restriction is the Poisson pencil of the open Toda lattice. 
In practice, its matrix representation in the reduced coordinates is obtained by 
the matrix representation of the periodic Poisson pencil (|3.2|) deleting the n-th 
row and column, and setting a n = in the resulting matrix. The Lax matrix of 
the open Toda Lattice (as well as the Hamiltonian functions) is obtained simply 
by setting a n = in the Lax matrix ( p.6|) of the periodic problem. In particular, 
the single polynomial Casimir of the Poisson pencil of the open lattice is obtained 
as ff^ 1 ) = ifW| o =0 . The open Toda lattice is thus a complete GZ manifold of 
type 2n — 1. 
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4 Separation of variables 



In this Section we will show that the Toda lattice fits the scheme described in 



Subsection JA_. We will follow the path of Section |3|, considering at first the 
periodic lattice, and then stating the suitable changes to be done in the open 
case. 

The periodic Toda lattice is a GZ manifold of dimension 2n and type {1, 2n — 
1}. Thus the rank of the transversal distribution Z must be 2, and the dimen- 
sion of the reduced PN manifold 2n — 2. We divide the procedure outlined in 
Subsection |2j] in three steps. 

Step 1. The transversal vector fields Z\ and Z 2 - 

Proposition 4.1 The vector fields Z\ = d bn and Z 2 = d an /(a\ ■ ■ -a n _i) satisfy 

Lie Zl P x = Z x A Y 1A , Lie Z2 P x =Z X N Y 2jl , (4.1) 
with Y ltl = a n _ 1 <9 an _ 1 - a n d an and F 2 ,i = -<9 6l /(ai ■ ■ ■ a n _i). 

This property is proven making use of the standard formulas for the Lie derivative 
of a bivector. 

Step 2. The action of on the Casimirs and the Xj coordinates. 

To discuss this issue it is useful to recall the expression of the second Casimir 
#(2) = (-1)™+! ai • • • a n and to expand det(£) with respect to the last column to 
get: 

det(£) = fi(b n - X)L nyn + / u 2 L n , n _ 1 + (-l) n+1 a„L nil , (4.2) 

where Lij are the determinants of the suitable minors of C. Taking into account 
the specific form of these minors one can easily see that it holds 

Proposition 4.2 The second Lie derivatives of and with respect to Zi 
vanish. Furthermore, Liez^H^) = and Liez 2 (H^) = 1, so that the matrix 
G b a = Liez a (H ( - b \\)) introduced in Subsection f£J| has the form 



G(X) 



( Lie Zx {HW) 0^ 
v Lie^HW) 1 / 



(4.3) 



Thus the Xj coordinates are the roots of the monic degree (n — 1) polynomial 
Lie Zl (H {1) ) = d bn det(C) =L n>n . 

Step 3. The action of the vector field Y and the fij coordinates. 

We have to consider the vector field Y = Y\ t \ = a n -i<9a n _i — a n d an , and to discuss 
its action on the exact eigenvectors of N*. According to the discussion following 
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Proposition [27|, to construct such an eigenvector we must find a vector (l,p(A)) 
such that (1, p(A))G(A) = for A = \j. Since this vector is simply given by (1, 0), 
we have that H^ 1 ' is an exact eigenvector of N*, and this is true, for all r, for 
Y r (H^) as well. In order to build a normalized exact eigenvector, we have to 
analyze a bit further the terms in the expansion ( |4.2j ) of the determinant of C 
Actually, one has that: 

1. L n n is independent of a n and a n _i; 

2. a„L ni i = +CK± where K\ is linear in a n and does not depend on a„_i 
and /z; 

3. fi 2 L n n _i = C 2 + CK 2 where K 2 is linear in a n _i and does not depend on 
a n and fi. 

Thanks to the linearity properties of Kj, we have that Y(H^) = Y(Ki + K 2 ) 
satisfies the recursion property Y 3 (H^) = Y(H^). This ensures that the func- 
tion 

F = log(Y(#«) + Y 2 (H^)) 

satisfies Y(F)=1, and, according to Theorem [2.4| , is the desired generator of 
the fj,j coordinates. We notice that, due to the cyclic nature of the periodic 
Toda system, the pair Z\, Z 2 of deformation vector fields is by no means unique; 
other admissible pairs can be obtained via a cyclic permutation of the indices 
(with a corresponding change in the vector field Y). It would be interesting to 
compare the DN coordinates considered here with the separation variables for 
Toda systems used in, e.g., [|, |K| |15|| . 

Finally, we state the corresponding results for the open Toda lattice. We have 
to look for a single deformation vector field, say Z\ it is still given by 9& n ; the 
vector field Y now is given by a n _i9 an _ 1 ; the recursion relation on Y(H^) closes 
at the first step, Y 2 (H^) = Y(H^), and the generating function can be taken 
as F = log (Y(H^)y 



5 Example: the three— particle case 



We close the paper with some explicit expressions for the three-particle case. 
The Poisson pencil is written in equation ( |3.2| ). The Lax matrix is given by 



/ 



C 



\i (&i — A) — fx 2 
a 1 [i (b 2 - A) 



«3 



\ 



a 2 



li (h - A) 



(5.1; 
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The spectral curve is 



with 



-C 2 - (A 3 + H \ 2 + H x \ + H 2 )C + K = 0, 

H = -(h + b 3 + b 2 ), H x = a 2 + b x b 3 + a x + b 2 b 3 + a 3 + b x b 2 , 
H 2 = -(bib 2 b 3 + a x b 3 + b x a 2 + a 3 b 2 ), K = a x a 2 a 3 , C = /z 3 . 



There are two nontrivial flows, given by: 

di = a\b 2 — a\b\ a x = aibib 3 + a x a 3 — aib 2 b 3 — a\a 2 

Xi — \ b\ = a 3 — cl\ X 2 = < bi = a\b 3 — 0362 

and cyclic permutations and cyclic permutations 

(5.2) 

The transversal vector fields are Z x = db 3 and Z 2 = d a3 /a x a 2 , and we have 
Y = a 2 d a2 — 03^3. The DN coordinates can be found as follows. The roots of 
the polynomial 

Lie Zl (H (1 \X)) = d b3 (H Q X 2 + H X X + H 2 ) = -A 2 + (61 + b 2 )\ - (b x b 2 + a x ) 

are Ai and A2. Then [i x and [i 2 are given by the function 

F = log (Y(H^) + Y 2 (H^)) = log(2a 2 A - 2a 2 b x ) 

evaluated at A = Ai, A 2 . 

For the open case, the Poisson pencil can be computed from Q3.2|), according 
to the procedure outlined at the end of Section [| 



popen 
X 







a x a 2 



—a x a 2 {pi — A) a x (A — b 2 ) a x 







(b 2 - A) a 2 (A - b 3 ) a 2 



(X-b x )a x 
(b 2 - A) a x (A - b 2 ) a 2 
(63 - A) a 2 



-a x 



a 2 



-a 2 



(5.3) 



and the spectral curve is the rational curve 

fi 3 = C = -A 3 + (b 2 + b x + 63) A 2 - (6i6 2 + a x + b 2 b 3 + 6163 + a 2 ) A 
+ b x b 2 b 3 + a x b 3 + b x a 2 . 

The separation variables can be constructed, mutatis mutandis, as in the periodic 



case. 
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